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would seem to be two alternative ways of approaching this
result (both briefly indicated in my previous article ([101],
p. 91). Their equivalence may be proved as follows.

Let w0 denote the summation operator [i, i, ..., i].
Then the simple unweighted sums of the columns of R
may be denoted by w0R, a row-vector like w0 itself.
Similarly, the sums of the columns of R* (say) may be
written (i) w0(R X R X R X R} ; but this is clearly the
same as (ii) (w0R) R X R X R.

The identity suggests two different working procedures;
(i) we may begin by multiplying the matrix R by itself :
repeat the process m times; and then end by adding the
columns of the final product-matrix ; or (ii) we may begin
by adding the columns of R, multiply the whole matrix
by this row-vector, i.e. weight each row of correlations and
add once again; then, if we repeat the multiplications m
times, the final sums will obviously be the same as before.
We may term (i) ' table-by-table multiplication' and
(ii) * table-by-column multiplication.'

Hotelling's iterative method might be described as a telescoped
form of (ii), and his matrix-squaring method as a curtailed form of
(i). By taking guessed approximations at each step instead of
exact multipliers, the final result is more rapidly reached than if the
full figures are retained at every step ; but at the same time the
underlying nature of the process is obscured. When the computer
carries out the calculations in full, he quickly perceives that the
differences between the successive approximations are each, the same
fraction of its predecessor, or nearly so : hence there is no need to

mathematicians for solving characteristic equations and the like. My late
colleague, R. C. Rowland, for example, has described it in his Note on a
Type of Determinantal Equation [62] in connexion with the solution of
certain dynamical problems in physics. I do not think, however, that it had
been explicitly shown before that the ultimate result would invariably be a
* hierarchy J: a proof without matrix notation was given in [102], p. 183.
-Since, by the product-moment formula, R = MM' (where M is the initial
matrix of measurements), the columns of R (which in its most general form
may consist of covariances, averaged or unaveraged, adjusted or unadjusted)
may be described as comprising the * first moments,' and the columns of
R* = MM'MM' . . . MM' as containing * higher moments' (cf. [93],
p. 287; [102], p. 178). Hence these terms have occasionally been used to
distinguish the two methods.